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Abstract. We determine the two loop corrections to the Green's function of a quark current 
inserted in a quark 2-point function at the symmetric subtraction point. The amplitudes for the 
scalar, vector and tensor currents are presented in both the MS and RI'/SMOM renormalization 
schemes. The RI'/SMOM scheme two loop renormalization for the scalar and tensor cases 
agree with previous work. The vector current renormalization requires special treatment as it 
must be consistent with the Slavnov- Taylor identity which we demonstrate. We also discuss 
the possibility of an alternative definition of the RI'/SMOM scheme in the case of the tensor 
current. 
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1 Introduction. 



Non-abelian quantum field theory underlies the strong nuclear force which binds quarks and 
gluons into hadrons. At high energy these quarks and gluons are asymptotically free, [H [2], 
and so a good approximation to the physics of hadrons in deep inelastic scattering can be 
achieved by perturbation theory. When the coupling constant, g, is small then the only difficulty 
is the actual computation of a large number of Feynman diagrams which prevents one from 
obtaining precise estimates. However, the physics of the specific structure of hadrons resides 
in the non-perturbative or low energy regime where, because the coupling constant is large, 
then perturbation theory is not applicable. Instead one focuses on the computation of matrix 
elements involving the relevant operators for the hadrons or deep inelastic scattering process. 
In principle such matrix elements can be measured accurately by using a lattice regularization 
of the non-abelian gauge theory. If one has access to powerful enough computers then one can 
build a solid picture of the dependence of the matrix elements with momentum scales. One issue 
which arises in the lattice computations is that whilst concentrating on the low energy aspect, 
the resulting matrix elements must still match onto the high energy behaviour which one can 
calculate in perturbation theory. This is not a trivial exercise. For instance, the operators one 
has to consider undergo renormalization. In perturbation theory the anomalous dimensions of 
key operators are known to at least three loops. See, for instance, [3j HI O El [8j [9] . However, 
this is invariably in the standard (non-physical) renormalization scheme known as MS. In 
this scheme essentially only the basic infinities with respect to the regularizing parameter are 
subtracted leaving the finite parts unsubtracted in the remaining part of the Green's function. 
Invariably one uses dimensional regularization in d = 4 — 2e dimensions. Therefore, the finite 
parts of the matrix elements at high energy are a reflection of the scheme. 

By contrast, the lattice uses various renormalization schemes which are different and physical. 
So to perform any matching in the overlap region requires knowledge of the matrix element in the 
same scheme whether this is MS or a lattice based scheme. The latter set of schemes are chosen 
primarily to reduce the financial cost of any numerical evaluation. For instance, derivatives 
within an operator or at any point require more computation. So the suite of lattice schemes 
are designed to minimize such complications. In earlier work the regularization invariant (RI) 
scheme and its modified version (RI') were defined in lattice computations, [10] 111] , and later 
developed to three and four loops for the continuum in several articles, both for the Landau 
gauge, [12], and general linear covariant gauges, [13] . Indeed matrix elements for deep inelastic 
scattering operators were evaluated to three loops in RI' in [13] [14] [T5] . More recently a variation 
on the RI' scheme has been developed, [16]. This is specifically related to matrix elements and 
designed to overcome a problem with potential infrared singularities. In essence the RI' scheme 
for 3-point and higher Green's functions involves subtracting the divergences at an exceptional 
momentum configuration. In other words the operator insertion is at zero momentum. To avoid 
this exceptional point and hence the related infrared issue, the RI'/SMOM renormalization 
scheme was introduced, [16]. The annotation indicates that there is non nullification of any of 
the external momenta in a 3-point function. Indeed the external momenta are non-zero and 
their squares are all fixed to be at the same value when the Green's function is renormalized. 
Hence one refers to it as a symmetric subtraction point. 

Initially this scheme was applied to the scalar, vector and tensor currents at one loop, [16] . 
and to the scalar and tensor at two loops in |17] I18], More recently it has been extended to various 
low moment operators used in deep inelastic scattering at one loop, [19]. However, the main 
focus of [17] [T8] was the construction of the anomalous dimensions and thence the conversion 
functions from the RI'/SMOM scheme to the MS one. These are central to any mapping of 
lattice results to the perturbative region for measurement comparisons. However, when one 
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undertakes any measurement on the lattice the Green's function with the operator insertion 
has free Lorentz indices and therefore it has to be written in terms of a set of basis tensors. 
Then different combinations of the free Lorentz indices can be determined and information 
on the associated scalar amplitudes extracted. Whilst \17\ [T8] concentrated on the operator 
renormalization, the explicit values of the individual two loop scalar amplitudes were not given 
which would be invaluable for lattice measurements. Therefore, it is the main purpose of this 
article to provide that information at two loops not only for the Green's function with scalar and 
tensor current insertions to augment the work of [171 [T8] but also for the vector current. This 
is partly because the latter has not been treated within the RI'/SMOM formalism at two loops 
but mainly because it underlies the renormalization of the deep inelastic scattering operators 
considered at one loop in [19] . The amplitudes for those two sets of operators will be considered 
separately, [20]. The case of the vector current is special as its renormalization is connected 
with the Slavnov- Taylor identity as discussed in [16]. However, one feature of the tensor current 
is that the actual definition of the RI'/SMOM scheme for such operators is not unique. This is 
because there is a relatively large set of basis tensors due to the free Lorentz indices. Therefore, 
a different basis choice would lead to different scheme definitions as we will indicate. In addition 
the way one projects out the part of the Green's function whose finite part will be absorbed 
into the operator renormalization constant is also subject to a large degree of choice. As there 
is a range of ways of defining the RI'/SMOM scheme for the tensor current we will give one 
alternative for illustration but will also present all the amplitudes for the currents considered 
here in the MS scheme too. So an interested reader has the liberty to toy with variations on the 
RI'/SMOM scheme definition of [TBI \T7\ [T8] . 

The article is organized as follows. General aspects of the computations we perform as 
well as the techniques used to carry them out are given in section two. The three specific 
operators we consider and the details associated with each are discussed in the three subsequent 
sections. Aspects of the conversion functions are given in section six including that for the 
alternative scheme devised for the tensor current to allow one to contrast with that of [17] [15] . 
Our conclusions are given in section seven. The main results are presented in a series of Tables. 

2 Preliminaries. 

To start with we focus on the generalities of the computation we are interested in. The three 
basic quark currents are the scalar, vector and tensor currents and we use the compact notation 
introduced in [T91I2T] 

S = tpip , V = rp^Tp , T = rpcr^ip (2.1) 

where a^ v = |[7^, J v ] is antisymmetric. Each of these operators, O, is inserted into a quark 
2-point Green's function where the two independent external momenta, p and q, flow in through 
each external quark leg as illustrated in Figure 1. As we will be concentrating on the renormal- 
ization of the operators and the consequent finite parts of the Green's function of Figure 1 at a 
symmetric renormalization point, we note that from here on we take, |16l 117} H8], 

P 2 = q 2 = {p + q? = - V 2 (2-2) 

which implies 

pq = \i? (2.3) 

where \x is the mass scale associated with the renormalization point. At this point the Green's 
function can be written in terms of a set of scalar amplitudes with respect to some basis of 
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Lorentz tensors. Even though the scalar current has no free Lorentz indices there are two 
independent amplitudes as the two independent momenta lead to two independent structures 
deriving from the 7-matrices. This was discussed in \19\ but we note that when either of the 
external momenta is nullified, equating to the RI' momentum configuration, then there is only 
one independent tensor in the basis. Therefore, in general we can write the Green's function for 
each of our operators 0* ^ separately as 

rii 

k=l 

where i is the label corresponding to the operator, which is either S, V or T here. The scalar 
amplitudes are denoted by q) with Vn^^ ^ {p,q) corresponding to the basis tensors. 

The latter were given in [19j and were derived by writing all the one loop Feynman diagrams 
in terms of a basic set of master tensor integrals. In other words each diagram was stripped 
of all features to leave purely basic integrals. These were evaluated by standard methods and 
then substituted back so that the computation could be completed by contracting the Lorentz 
indices with the stripped off 7-matrices. The reason for proceeding in this way was to ensure 
that no tensors contributing to the basis in the decompostion were omitted. Ordinarily one 
would construct the basis by building all possible tensors from the basic tensors such as tw, 
7-matrices and independent momenta. Although here there are two momenta it turns out that 
even for simple currents the basis can be quite large. This is aside from the fact that we are 
only constructing the basis at the symmetric subtraction point where the values of p 2 , q 2 and 
(p + q) 2 are au equal. Away from this point the basis of tensors will be significantly larger but 
we are not interested in the form of the amplitudes not at the symmetric point. 



t p+q 




Figure 1: Momentum flow for the Green's function, ( ip(p)0' t tll ^ n (— p — q)ip(q)). 



As the scalar amplitudes are the quantities we seek then it is possible to compute each 
individually via a projection onto the Green's function itself. In other words 



= MUvtrT^M ({mo^..., n x-p - qmq) 



(2.5) 



where M\i is a matrix of rational polynomials in d with k and I labelling the different amplitudes 
relative to the ordering of the tensors in the basis, |19j . Thoughout we use dimensional regular- 
ization in d = 4 — 2e dimensions. This matrix is computed from the actual tensors themselves. 
First, if we define the matrix 



Ml 



ki 



^lk)^ n Xp^)^\i)^ nr {p^) 



(2.6) 
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then is its inverse. Of course the choice of basis tensors is not unique and we use the one 
constructed for each operator at one loop. Whilst that involved stripping the spinor structure off 
all the Feynman graphs and then substituting the master tensor integrals, such an exercise would 
be too difficult and cumbersome to implement at two loops. Therefore, we use the projection of 
the Green's function as outlined here. This was also used at one loop and reproduced the direct 
evaluation results of |19j . Therefore, we are confident that we have a complete basis. Though 
there are additional checks on the results such as agreement with the Slavnov- Taylor identity, 
in the case of the vector, as will be discussed later. 

With our general decomposition of the Green's function we can now discuss the renormal- 
ization and the method to define the RI'/SMOM scheme renormalization constants. Of the set 
of scalar amplitudes, one or more will contain the poles in e. If we denote this amplitude, or set 
of amplitudes, by the label then the renormalization constant for the operator is defined by 
the condition, [TBlfTTlfTS]. 



lim 



„W' 7 M'/SMOM o , 



1 . 



(2.7) 



In other words there are no 0(a) corrections after renormalization where we set a = g 2 /(167r 2 ) 
and g is the coupling constant. In addition we note that the origin of the second aspect of 
this scheme definition is that the quark wave function renormalization is carried out in the 
RI' scheme. Briefly, this scheme is defined by ensuring that after renormalization there are no 
0(a) corrections in 2-point functions but for 3-point functions and higher the renormalization 
is performed as one does in the MS scheme. This RI' scheme and its sister scheme, RI, were 
introduced in |101 lll| and examined in the continuum case at three and four loops, [12 \ I13 |. ITS] . 
The latter work was for a linear covariant gauge fixing. Such computations are important when 
one is converting from schemes such as RI' or RI'/SMOM to MS as one has to express the 
parameters of each scheme in terms of the same parameters of the other scheme. It transpires 
that the coupling constants are in direct correspondence, [15] . 



2 MS + °( ( 



MS) 



(2.1 



to the order they were computed where we indicate the scheme to which the variables relate via 
the subscript. However, for a linear covariant gauge the gauge parameter, a, is not in one-to-one 
correspondence except in the Landau gauge as, [15] . 



a 



RI' 



1 + 



-9a- 



MS 



18q MS - 97) CU + S0T F N f 



a- 



MS 



•MS - j ■ ■ ■ -v 36 
+ ((l8o^g - 18a|jg + 190a|jg - 576^3)^ + 463^ + 864C(3) - 7143) C\ 



+ 



+ 



(- 320a|jg - 320 aM g + 2304C(3) + 4248) C A T F N f 

a m + °( 



4608C(3) + 5280) C F T F N f ) 



a Ms) 



(2.9) 



to two loops. The full three loop result is given in [13J. These relations between the parameters 
of each scheme are required when converting the amplitudes. The group Casimirs are defined 
in the usual way by 



rpCLrpa 



C F 



acd fbcd 



f aca f 



C A S 



ab 



Tr ( T a T 



T F 5' 



ab 



(2.10) 



where T a are the generators of the colour group whose structure functions are f abc . Throughout 
we work with Nf massless quarks so that all our expressions for amplitudes are for the chiral 
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limit. In principle quark masses could be included. However, the basic master scalar Feynman 
integrals at two loops with massive propagators are not known exactly for the symmetric point. 
So quark mass dependence for the scalar amplitudes of the Green's function we are computing 
could only be determined in, say, the small quark mass limit. In connection with the relation 
of the parameters between schemes we should note that in [16] a different parameter mapping 
was chosen. There it was assumed that the gauge parameter was the same in both schemes 
in the same way that the coupling constants are. Here we will present all our results for the 
non-Landau schemes in what we regard as the full RI' context which led to (|2.9p . 

The next aspect of (|2,7p which we need to draw attention to is the way in which the operator 
renormalization constant is actually defined. For instance, in writing (12. 7p we are making a 
basis dependent statement. The choice of the basis tensors is purely arbitrary and another 
choice with the same underlying criterion of having no 0(a) correction will clearly lead to a 
different numerical value of the corrections in the renormalization constant itself. Moreover, 
the amplitudes will also have the same degree of arbitrariness. This is not the same as scheme 
dependence. In that situation if one could compute the amplitudes to all orders the physics 
would not be affected by the renormalization scheme choice. However, as one has to truncate 
series in quantum field theory due to the limit of calculability it is clear that the construction of 
conversion functions, such as those of [16j [T71 [18] , could be affected by the basis choice. Though 
in [THJ, [171 [18] the approach used was to project the Green's function with the appropriate Born 
term before rendering that projection to have no 0(a) part. This defined the RI'/SMOM scheme 
for the scalar and tensor currents. However, once the full scalar amplitudes have been computed 
it is clear that there is in fact a sizeable number of ways of defining a so called RI'/SMOM scheme. 
This was discussed in the one loop context in [19] where it is becomes a more important issue for 
the operators used in deep inelastic scattering. Therefore, we will discuss a possible variation on 
the scheme of [T71 [18] for the tensor current at two loops. In addition, partly because of these 
general aspects we will provide all the amplitudes in the reference MS scheme. This is primarily 
for lattice computations where it is easier to convert to MS on the lattice before looking at the 
matching onto the ultraviolet part of the Green's function where perturbation theory is valid. 

Next we address some of the technical aspects of the computation. One feature of the basis 
of tensors is that they involve products of 7-matrices. This is because one can have contractions 
such as j$ and Therefore, we have chosen to work with the generalized T-matrices of [22l t23[ [24] 
which are denoted by T^"^" . They are defined by 

pMl...Mn = 7 [mi _ _ _ 7 Mn] (2.H) 

which is totally antisymmetric in all the Lorentz indices and the notation includes the overall 
factor of 1/n!. So, for instance, = General properties are given in |25[ 126] . One 

advantage of this choice for the tensor basis is that in T-space 

tr (rf^'T^p) oc 5 mn I^-^- M " (2.12) 

where l^---^m,v\...v n j g ^ e unit matrix in this space. So there is a natural partition for the basis. 
The main tool to handle the tedious algebra for manipulating tensor projections is the symbolic 
manipulation language Form, [27] , The Feynman diagrams are generated with the Qgraf 
package, [28J, and converted into Form notation whereby the colour and Lorentz indices are 
appended. For the Green's function we are interested in with quark current insertions there are 
1 one loop graph and 13 two loop graphs to be computed. These graphs are readily broken up 
into a set of tensor integrals where the external momenta swamp any free Lorentz index when 
we have multiplied by the appropriate tensor of the projection basis. As these integrals are in 
essence 3-point functions evaluated at the symmetric point the next stage is to break them down 
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into the known scalar master integrals. These are collected in [16] but were derived in various 
articles, [29j [301 Ell E2]- The engine room of this aspect of the computation is the Laporta 
algorithm, [33]. Briefly, that method allows one to build a redundant set of equations where 
basic Feynman graphs with irreducible numerators are all related by integration by parts and 
Lorentz identities. These can then be solved as a tower of equations with the master integrals 
being the foundation for the many sets of integrals in predefined sectors. There are a variety 
of computer packages available to build a Laporta system, [34], I35j . However, we have used 
Reduze, |36j, which uses Ginac, [37J, and involves C++ at its root. For the particular Green's 
function we consider we need only build two topologies using the Reduze package since all 
the Feynman graphs are either based on a two loop ladder or the non-planar ladder. Once the 
system is built it can be converted into Form language and all the two loop tensor integrals for 
all the Feynman graphs at the symmetric subtraction point are written in terms of the known 
scalar master integrals listed in [18] . A check on the Reduze results is that we do reproduce the 
expressions given in [16], H7J H8] . Though the expressions we give at two loops for the amplitudes 
are new. 

As was evident in |17t 118] the final forms for the two loop anomalous dimensions and as- 
sociated conversion functions were surprisingly long. As we will be presenting a large number 
of amplitudes for the three currents in two renormalization schemes in order to save space we 
will collect the main results in Tablefl To do this we have had to split the amplitudes by their 
colour group structure and so we have defined 



where the parameters are in either scheme. Here we will adapt the convention that all expressions 
are in the RI'/SMOM scheme unless explicitly indicated to be in the MS scheme. The coupling 
constants are the same but the gauge parameters are not, (12. 9h . In (I2.13j) the entities al denote 
the basis of numbers which appear in that specific part of anomalous dimension. This includes, 
for instance, the gauge parameter dependence and it is evident from the left hand column of 
each table what these numbers actually are for each colour structure. In essence they relate to 
the form of the scalar master integrals at the symmetric subtraction point. The other entities, 

c (i)n ' are ^ ne actual coefficients of the number basis. The summation label n corresponds 
to the row of each table. The superscript k annotates the loop order, as the first number, 
and at two loop the second number references the colour group structure. Although the exact 
expressions at two loop represent all the information on the scalar amplitudes which we seek, 
we have relegated the Tables to the end of the article as it is the numerical evaluation which is 
ultimately of practical use. These expressions will be provided for each operator in succession 
in the next few sections. Throughout the one loop amplitudes are in agreement with those of 

EE]- 

3 Scalar current. 

In this section we concentrate on the scalar current. The aim for this and the other currents 
is to provide not only the anomalous dimensions in the full RI'/SMOM scheme but also the 

* Attached to this article is an electronic file where all the expressions presented in the Tables are available in 
a useable format. 
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finite parts of all the amplitudes. Previously in |17} [T8] the two loop anomalous dimension 
and conversion function were given. However, in order to assist with the extraction of results 
from the lattice, measurements have to be made in different components of the Lorentz basis of 
tensors and therefore all the finite parts of the Green's functions need to be known accurately. 
As the amplitudes for the scalar current have not been given at two loops we briefly report on 
this situation in this section first. The decomposition into the projection tensors involves two 
tensors which are 



(0) 



P(2)(P>9) 



1 



2 (2) 



(3.1) 



• ■(Un 



where we use the convention that if a momentum is contracted with a Lorentz index of T?\ 

(n) 

then we replace that index by the momentum to save space. The matrix used for constructing 
the explicit projection is given by, [19], 








(3.2) 



This is a simple example of the partitioning of the matrix as a consequence of the choice of the 
rV n ) basis. As noted previously to record the full two loop expressions for the amplitudes would 
be demanding on space. However, as there are only two amplitudes for the scalar case we present 
the results for both schemes in Tables 1 to 4. The notation of fl2. 13f> is used with the convention 
that the MS scheme results are annotated explicitly. Otherwise the results are in the RI'/SMOM 
scheme. For the scalar current the RI'/SMOM scheme renormalization condition is to project 
the Green's function with the Born term and then define the operator renormalization constant 
so that there is no 0(a) finite part, [TBI El OS]- Since we are using the generalized T-matrix 
basis then for the scalar case this equates to ensuring that there are no 0(a) corrections to the 
channel 1 amplitude. This is because of (|2.12p . Consequently there are no RI'/SMOM scheme 
coefficients for channel 1 in these Tables. 

Throughout we use the standard basis for the numbers which appear in this symmetric 
subtraction point momentum configuration for this Green's function, [TTj. There ip(z) is the 
derivative of the logarithm of the Euler Gamma function, 



Sn(z) = i=9 




(3.3) 



where Li n (z) is the polylogarithm function, £(z) is the Riemann zeta function and S is given by 
a combination of various harmonic polylogarithms, \17\ I32j. 



nfi + 'H 



(2) 
43 



(3.4) 



Whilst the exact two loop expressions are the output from the Form computation, for practical 
purposes expressing the result in numerical form will be more pragmatic for users. Therefore, 
we record this in an explicit equation for the case of SU (3) not only for the scalar operator but 
also for the other cases we consider later. For the MS scheme we have 



£?i)(P,9) 



s f 2 )(P^) 



MS 



MS 



-1- [1.1040618a + 0.6455188] a 

- 48.4885881 + 7.58534654a + 3.67844381a 2 - 6.3468728^/1 a 2 + 0(a 3 ) 
[1.0417366a - 1.0417366] a 

- 11.1668053 - 7.37922016a - 3.5592665a 2 - 0.4629940^/1 a 2 

+ 0(a 3 ) . (3.5) 
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Those for the RI'/SMOM scheme are 



£fi)(p,?) = - 1 + 0(a 3 ) 

T,f 2) (p, q) = [1.0417366a - 1.0417366] a 

- 10.4943447 - 16.2776049a - 3.97172981a 

+ [1.1574851a - 0.4629940] N f ] a 

Throughout we use the numerical values 



0.7813024a 3 



2 + 0(a 3 ) 



(3.6) 



C(3) 



S3 



1.20205690 

.1838167 , s 2 



£ = 6.34517334 

7T N 



10.09559713 



0.32948320 , s 3 



0.32225882 , s 2 



0.19259341 . 



7T 



0.22459602 



(3.7) 



Clearly there is a weak dependence on the number of quarks in channel 2. To reinforce an 
earlier point the gauge parameter a in these sets of expressions is the variable in the RI'/SMOM 
scheme. Though the coupling constant, a, is in the same scheme it is an exact map of the MS 
variable. Moreover, as with other RI'/SMOM scheme anomalous dimensions of gauge invariant 
operators the higher loop corrections will depend on the gauge parameter. This is because the 
scheme is a mass dependent one and not a mass independent one like MS. Finally, we note that 
using (|2.9p the RI'/SMOM scheme anomalous dimension is 



7 5 (a,a) 



RI'/SMOM 



3Ci?a 



+ 



(3a 2 + 9a + 66)V> 



1 1 n 

3' 



(2a 2 + 6a + 44)tt 2 



+ 



9a^ - 27a - 555 

2 



16ir z + 156 - 24V> (f) T F N f 



C F a 2 
18 



+ 0(a 3 ) (3.8) 



which agrees with [161 El HE] when a = where we will annotate the anomalous dimensions 
with the scheme. The a dependent terms differ because of the different ways the renormalization 
of a is performed. We have chosen to define the a renormalization using the RI' scheme of |13j 
whereas in [T7] the a renormalization is taken to be in the MS scheme. As ultimately lattice 
computations are in the Landau gauge this difference in definitions would only be important in 
non-Landau linear covariant gauges. 



4 Vector current. 

As the vector operator has not received attention at two loops, we devote this section to it 
in detail using the general notation discussed previously. First, the basis of tensors used to 
decompose the Green's function at the symmetric subtraction point is, 1191 . 



rfaM = q 4 . p&„(p.9) = q 4 - t%»M = i r (3)^ 5 (4-1) 

H 1 f* h 1 
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where the matrix used to perform the projection into the various amplitudes is, [19J, 





f 9 


12 


6 




6 




12 


o \ 




12 


16(d- 1) 


8(d- 


1) 


8(d- 


1) 


4(d + 2) 





1 


6 


8(d- 1) 


4(4d- 


-7) 


4(d- 


1) 


8(d-l) 





36(d - 2) 


6 


8(d- 1) 


A{d- 


1) 


4(4d- 


-7) 


8(d-l) 







12 


4(d + 2) 


8(d- 


1) 


8(d- 


1) 


16(d- 1) 







V o 


















-i2 ; 



(4.2) 



With this setup we have applied the computational algorithm to extract the RI'/SMOM renor- 
malization constant. To do this for the RI'/SMOM scheme requires a different approach com- 
pared to the other quark currents. The reason for this is that the vector current is a physi- 
cal operator and therefore its renormalization is already determined by general considerations. 
Specifically as it is physical its anomalous dimension is zero which is widely known in the MS 
scheme. However, once it is accepted that the renormalization constant is unity in one scheme 
then it is unity in all other schemes, [38] . Underlying this is the Slavnov- Taylor identity which 
relates the renormalization of the Green's function with the divergence of the operator to the 
renormalization of the quark 2-point functions. Indeed this was discussed in [13] for the RI' 
scheme and demonstrated to be consistent to three loops. The situation for the RI'/SMOM 
computation is more involved as there are two momenta flowing through the Green's function 
with the operator insertion. Therefore, to reproduce the Slavnov- Taylor identity the Green's 
function of Figure 1 has to be contracted with the vector (p + This is one reason why we 
have to decompose the Green's function into a basis of projection tensors. Once we have estab- 
lished this then the contraction can proceed. For the case of the RI' scheme this aspect of the 
reconciliation with the Slavnov- Taylor identity is simplified significantly because there is only 
one momentum flowing through the Green's function. Indeed put another way the contraction 
of the graph of Figure 1 with (p + q)^ will effectively become the renormalization condition for 
Z v but will naturally produce unity as expected from general theorems. Whilst we are focusing 
on this feature for the vector current in detail here it transpires that the same issue arises for 
operators which are used in deep inelastic scattering. This was discussed at one loop in [19] for 
the sets of operators labelled W2 and W3 in the notation of [21]. This is because the contraction 
of Lorentz indices of one of the operators in each set is equivalent to the divergence of the vector 
current. Therefore, the renormalization of those operators also has to be consistent with the 
Slavnov- Taylor identity. 

To focus on the issue we concentrate on the MS renormalization first. The explicit results 
for each of the amplitudes is given in Tables 5 to 8. The numerical values for the amplitudes for 
SU(3) are 



s (i)(P. 



S (6)(^«) 



= - 1 + [1.6249301 - 0.5831936a] a 



MS 



MS 



MS 



+ 



6.1248321 - 3.2229010a - 1.8119992a 2 + 0.2362586iY, 



a 2 + 0(a 3 ) 



[1.4720824 + 0.3056953a] a 



+ 



18.7974908 + 4.5957818a + 1.0954082a 2 - 1.3299518-/V/ a 2 + 0(a 3 ) 



[1.7777778 + 1.1945842a] a 



+ 



MS 



44.3805855 + 12.1090504a + 4.2805934a 2 - 2.8641975^/ a 2 + 0(a 3 
2.0834731a 
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- 39.7873696 + 0.3484662a + 0.1736228a 3 
+ 0(a 3 ) . 



3.009461 IN, 



(4.3) 



As noted at one loop in [19] there is a degree of symmetry which is also evident at two loops. 
This is primarily due to the way in which we chose our basis of Lorentz tensors which is of course 
not unique. The fact that various channels pair off to two loops can be regarded as an internal 
check on both the construction of the projection matrix and the symbolic manipulation code 
used to perform the computation. In the tables for the vector case we do not reproduce columns 
for the channels 4 and 5 because of the above equalities which hold exactly to two loops. In 
order to see that the Slavnov- Taylor identity is satisfied to two loops in the MS scheme we have 
computed that combination of the amplitudes which correspond to the Green's function of the 
insertion of the divergence of the vector current. With the contraction with {p + q)^ there will 
be two terms. One will involve j$ and the other ^. For the former the combination gives 



1 



MS 2 (2) 
— 1 — ctC fcl + 

+ 0(a 3 ) 



s (2)(p>°) 



1 



MS 



MS 
13 



, x 41 , , 13 9 2 
3C(3)- T + 3C(3)a- y a--a 2 



C A + \c F + 7 -T F N f 
o A 



Cpa 1 
(4.4) 



and that for the latter is a different combination of amplitudes but produces the same result. The 
right hand side of ()4.4p is clearly the finite part of the quark 2-point function after renormalizing 
in the MS scheme. Therefore, given the fact that the Green's function is symmetric under 
interchange of p and q then this shows that the Slavnov- Taylor identity is correctly embedded 
within our computation with the vector current having a renormalization constant of unity. 

The situation for the RI'/SMOM scheme amplitudes is completely parallel to that of MS. 
We have given the amplitudes for this case in Tables 9 to 12 where we have used 

Z v = 1 + 0(a 3 ) (4.5) 

for the renormalization of the vector current. As a brief summary the 577(3) numerical values 
are 



1 + [1.6249301 + 0.7501398a] a 



+ 



31.5890381 + 12.2494724a + 2.8130241a 2 + 0.5626048a 3 
- [2.0970747 + 0.8334886a] N f ] a 2 + 0(a 3 ) 



S( 5 )(Pi9) 

[1.4720824 + 0.3056953a] a 



+ 



18.7974908 + 5.1040424a + 1.1463575a 2 + 0.2292715a 3 



[1.3299518 + 0.3396615a] N f ] a 2 + 0(a 3 



[1.7777778 



1.1945842a] a 



+ 



44.3805855 + 19.3949024a + 4.4796908a 2 + 0.8959382a 3 



- [2.8641975 + 1.3273158a] Nf] a 2 + 0(a 3 
2.0834731a 



- 39.7873696 
+ 0(a 3 ) . 



2.4294979a + 0.1736228a 2 - 3.0094611iY r 



(4.6) 
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In order to see that the definition of Z^ 1 is consistent with the Slavnov- Taylor identity we have 
repeated the computation of (|4.4p for the RI'/SMOM amplitudes. In this case the piece involving 
$ produces 

Eftfo?) - \^\ 2) {p,q) ~ = ~ 1 + 0(a 3 ) . (4.7) 

Here there are no corrections which is consistent with the RI'/SMOM scheme since, [16], it 
uses the RI' quark wave function renormalization. This is chosen in such a way that the quark 
2-point function has no 0{a) corrections. In other words the finite part of that 2-point function 
is absorbed into the finite part of the wave function renormalization constant. Therefore, the 
vector current anomalous dimension 



7 y (a,a) 



= <3(a 3 ) (4.8) 

RI'/SMOM 



is consistent with the underlying Slavnov- Taylor identity. Whilst 7 (a, a) , is zero to 

J ' v ' ' RI'/SMOM 

all orders we include the order symbol merely to record the order to which we have explicitly 
verified this to. 



5 Tensor current. 



Next we record the parallel situation for the tensor operator noting only the points where there 
are differences from earlier work and discussions. The Lorentz tensor basis for the amplitude 
decomposition is, [19] , 

1 



V, 



(p,q) 



(p, q) 



V, 



T 



1 



(3)/«/(P>9) - ~2 
1 



^(5)a«/(P»9) 



F(2) /zpPi/ 
r (2) liqPv 



r (2) i/</P/* 



[P/^ - Pf?A»] r (0) i 
^(4)^(P>9) = 72 

Ptow(p>?) 



(2)ti P qv 



(2) W ?f 



r(2)i/p9/i 

r (2) i^M 



(2) p 9 P/Ltft/ - 1 (2) pgPf #/Li 



P. 



(8)^1/ 



1 



2 ^( 4 ) » v p<i 



The matrix to determine the explicit projection is defined as, |19| . 



M 1 

with the four submatrices given by 

/ -9 



1 



M T 



11 



36(d-2)(d-3) 







Mb M 
M T 21 M T 



12 

T 
22 



\ 



M T 



12 



m\\ 



-12 -6 
6(d-2)(d-3) 

-12 -8(d-l) -4(d-l 

V -6 -4(d- 1) -4(2d- 5) / 

/ -6 -12 -12 -6 \ 



4(d-l) -2(d + 5) -8(d-l) 

\-2(d-\) -4(d-l) -4(d-l) / 

(-6 -4(d- 1) -2(d- 1) \ 

12 -2(d + 5) -4(d-l) 

12 -8(d-l) -4(d-l) 

\ / 



(5.1) 



(5.2) 



12 



M\ 2 



/ -4(2d-5) -4(d-l) -4(d-l) \ 

-4(d-l) -8(d-l) -8(d-l) 

-4(d-l) -8(d-l) -8(d-l)(d-2) 

\ 12 / 



(5.3) 



This produces the two loop MS numerical values for the amplitudes 



tfx)(p,q) 



s (8)(p>?) 



MS 



MS 



- 1 + [0.0623253 - 0.0623253a] a 

+ [l7.0099539 + 0.6409422a + 0.0544455a 2 - 1.6001145^1 a 2 + 0(a 3 ) 
[3.1252097 + 1.0417366a] a 



+ 



76.2022091 + 10.8541168a + 3.9065121a 2 - 5.5559283iV/ 



MS 



+ 0(a 3 ) 

S 5)(P.9) m 
[0.3056953a - 0.3056953] a 

- 1.8195802 - 3.2067998a - 1.0954082a 2 - 0.3396615^/1 a 2 + 0(a 3 ) 



MS 



MS 



MS 



MS 



[0.1528477a - 0.1528477] a 

- 0.9097901 - 1.6033999a - 0.5477041a 2 - 0.1698307^/1 a 2 + 0(a 3 ) 
[0.6113907a - 0.6113907] a 



3.6391605 - 6.4135995a - 2.1908165a 2 - 0.6793229iV/ 
1.0417366a + 3.1252097] a 



a 2 + 0(a 3 



- 76.2022091 + 10.8541168a + 3.9065121a 2 - 5.5559283iV/ 
+ 0(a 3 ) 



(5.4) 



with the full expressions for each amplitude in the MS scheme recorded in Tables 13 to 16. Unlike 
the vector case there is no Slavnov- Taylor identity to be satisfied by the renormalization constant 
of the tensor current. For the MS renormalization the renormalization constant is already known 
and we have used that here. However, as noted in [19] there is a variety of ways of defining 
the renormalization in the RI'/SMOM scheme case. In [161 [T71 [T8] the renormalization was 
performed by first projecting the Green's function with the tensor of the Born term. For the 
tensor current this is r"/L. The RI'/SMOM scheme renormalization is then defined so that there 
is no 0(a) correction in this projection. However, there is no reason to regard this as the unique 
way of defining the renormalization constant. Now that we have the complete decomposition 
into the tensor basis, one could define the renormalization so that instead the coefficient of 
has no 0(a) piece after renormalization. This tensor channel contains the poles in e which have 
to be removed. To us this is also a perfectly reasonable way to define the scheme. Though 
it depends of course on the other elements of the tensor projection basis which is not unique. 
Indeed in [19] this alternative way was studied and the one loop correction was found to be 
numerically smaller than that of \T7\ [T8] . However it was not clear if this would persist to 
next order. If not then it may be possible to improve the convergence by redefining the basis. 

First, we record that we have followed the original RI'/SMOM scheme definition of [16 1 117 } [T8] 
and reproduced precisely the full two loop renormalization constant and anomalous dimension 
of [16} [TT1 118j . This acts as a non-trivial check on the Reduze database of integrals we have 
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constructed by the Laporta algorithm. However, the full RI'/SMOM amplitudes have not been 
presented and we record the numerical values are 



s fi)0> q) 



s (2)(p> q) 



-1 + [0.1528477a -0.1528477] a 

- 0.9426788 - 2.9046959a - 0.7440869a 2 - 0.1146357a 3 

- [0.1698307 — 0.1698307a] Nf] a 2 + 0(a 3 ) 
[1.0417366a + 3.1252097] a 



+ 



76.8746697 + 18.8265135a + 5.2449633a 2 + 0.7813024a 3 
- [5.5559283 + 1.1574851a] Nf] a 2 + 0(a 3 ) 



(6)1 

[0.3056953a - 0.3056953] a 

- 1.8853576 - 5.8093917a - 1.4881739a 2 - 0.2292715a 3 
- [0.3396615 - 0.3396615a] N f ] a 2 + 0(a 3 ) 



2(4) (P>?) = Ej 5) (p,q) 



s f 7 )fe q) 



s f 8 )(^ q) 



(5)1 

[0.1528477a - 0.1528477] a 

- 0.9426788 - 2.9046959a - 0.7440869a 2 - 0.1146357a 3 

- [0.1698307 -0.1698307a] iV)] a 2 + 0(a 3 ) 
[0.6113907a - 0.6113907] a 

- 3.7707152 - 11.6187834a - 2.9763478a 2 - 0.4585430a 3 

- [0.6793229 — 0.6793229a] Nf] a 2 + 0(a 3 ) 

- [1.0417366a + 3.1252097] a 

76.8746697 + 18.8265135a + 5.2449634a 2 + 0.7813024a 3 

- [5.5559283 + 1.1574851a] Nf] a 2 + 0(a 3 ) . 



(5.5) 



The full explicit results are given in Tables 17 to 20. Whilst the symmetry derived from the 
interchange of the momenta p and q is evident numerically here, this reflects the actual symmetry 
in the exact expressions which we have checked explicitly. So we have omitted those columns in 
the tables corresponding to the relations given in (|5.5p . To two loops the associated anomalous 
dimension is 



7 T (a,a) 



RI'/SMOM 



Cfcl 



+ 



(9a 2 + 27a - 66)^'(|) - (6a 2 + 18a - 44)tt 2 



+ 



9a' 



24^ 



27a + 1035 



3/ 



16tt z - 252 



C A - 513(7 F 

C F a 2 



T F N f 



54 



+ 0(a 3 ) (5.6) 



where 7 T (a, 0) , agrees with I16LI171H51 . As with the scalar current the difference with 

' v ' 1 RI'/SMOM 6 1 — * — M — ! 

the a dependent terms between (|5.6p and that of |17j resides in the different ways the gauge 
parameter is renormalized. We again use that derived in the RI' scheme, |13j . There careful 
attention was given to ensuring that the full three loop renormalization of QCD was consistent 
with the Slavnov- Taylor identities in that scheme. 
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We close this section by discussing an alternative way of defining the tensor current renor- 
malization which was noted in [19] . Instead of the approach of [T7] the finite part of associated 
with in the tensor basis was used to define the renormalization constant, with respect to 
the basis we have introduced. This channel contains the singularities in e which must always 
be absorbed in any scheme. Consequently we find the two loop anomalous dimension in this 
alternative scheme is 



7 (a, a) 



alt RI'/SMOM 



C F a 



+ 



(45a 2 + 135a - 330)V>'(§) - (30a 2 + 90a - 220)vr 2 



81aT - 243a + 3501 



C A 



120i//(§) - 80vH - 900 T F N f 



1539C F 
C F a 2 



162 



+ 0(a 3 ) . (5.7) 



For SU(3) the numerical equivalent is 



7 T (a,a) 



alt RI'/SMOM 



1.3333333a 



+ 



1.9065121a 2 + 5.7195362a + 40.9078004 - 1.9674761iV) 



+ 0(a 3 



(5.8) 



Clearly the leading term is the same as the original RI'/SMOM scheme as it ought to be. With 
this scheme choice the amplitudes are virtually the same as those for the original RI'/SMOM 
scheme of |19j which we have presented already. The only differences are that the channel 1 
coefficients in Tables 17 to 19 are absent whilst the coefficients of Table 20 are completely differ- 
ent. For this alternative scheme the appropriate coefficients are given in Table 21. Numerically 
we have 



s fl)0,9) 
s f 2 )(P^) 



alt 
alt 

alt 



- 1 + 0(a 3 ) 
[1.0417366a + 3.1252097]a 

+ [76.3969887 + 19.1449675a + 5.4041904a 2 + 0.7813024a 3 
- [5.5559283 + 1.1574851a]iY/]a 2 + <3(a 3 ) 



alt 



S (6)(P>9) 

[0.3056953a - 0.3056953]a 

- [1.8386328 - 5.7159421a - 1.5348987a 2 - 0.2292715a 3 
- [0.3396615 - 0. 339661 5a] N f ] a 2 + 0(a 3 ) 



sj)(p,g) 



alt 



alt 



S (8)(P>9) 



alt 



alt 



S (5)(P>9) 

[0.1528477a - 0.1528477]a 

- [0.9193164 - 2.8579710a - 0.7674493a 2 - 0.1146357a 3 

- [0.1698307 -0.1 6983070]^/] a 2 + 0(a 3 ) 

[0.6113907a - 0.6113907]a 

- [3.6772656 - 11.4318841a - 3.0697974a 2 - 0.4585430a 3 

- [0.6793229 - 0. 6793229a] N f ] a 2 + <3(a 3 ) 

- [1.0417366a + 3.1252097]a 
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- [76.3969887 + 19.1449675a + 5.4041904a 2 + 0.7813024a 3 

- [5.5559283 + 1. 157485 la] N f ] a 2 + 0(a 3 ) . (5.9) 



For the Nf independent part at two loops the numerical differences between these amplitudes 
and those of (|5.5p is insignificant. 



6 Conversion functions. 



In this section we record the conversion functions for changing from one scheme to another for 
the various operators considered here. For an excellent background to this see, for example, 
|38j . The conversion functions, Cj(a, a), are defined from the explicit forms of the appropriate 
renormalization constants in each scheme. So, for instance, 



C s (a,a) 



7 s 
RiySMOM 

z m 



C v (a,a) 



7 V 
RiySMOM 

Z MS 



C T (a,a) 



^RI'/SMOM 

z m 



(6.1) 



However, in deriving the explicit expressions in each case one must make a choice of scheme 
in which to express all the parameters themselves in. Here, we will use the MS scheme as the 
base for the variables. So that in the conversion functions, (|6.ip . a and a are MS parameters. 
Further, in practical terms this means that the RI'/SMOM scheme renormalization constant has 
been determined as a function of the RI'/SMOM scheme version of a and a. These have first 
to be mapped to their MS equivalent before the ratios in (|6.ip can be computed. Otherwise one 
might find that the conversion functions are not finite in four dimensions as they ought to be. 
With our definition of the RI'/SMOM scheme based on the RI' scheme definition of a, |13| . we 
have 



C 5 (a,a) 



1 



(3a + 9)V>'( 



(6 + 2a)vr 2 - 9a - 36 



Cfcl 



(72a 2 + 432a + 936)(t//(±)) 2 - (1248 + 576a + 96a 2 )^' a)n 2 



- (5616 + 864a + 432a 2 )V>'( ±) - (72 + 36aX'(§) 

- 15552s 2 (f ) + 31104s 2 (f ) + 25920s 3 (§) - 20736s 3 (§ ) 
+ (32a 2 + 288a + 608)vr 4 + (288a 2 + 576a + 3744)vr 2 
+ 648a 2 + 2592a + 1539 + (648a + 1944)S 



ln 2 (3)vr ln(3)vr vr 3 

2592C(3) - 108 y— + 1296 y ' + 116^= 

\/3 v3 v3 



C F 



192V/(i)vr 2 - 144(V/(i)) 2 + (162a 2 + 756a + 8226)^'(|; 



- (45 + 9a)ip'"(l) + 7776s 2 (f ) - 15552s 2 (f ) - 12960s 3 (f ) 
+ 10368s 3 (|) + (24a + 56)vr 4 - (108a 2 + 504a + 5484)vr 2 

- 486a 2 - 2268a - 34695 + (324a + 1620)S + 6480C(3) 



+ 54 



960vr< 



ln 2 (3h 
1440V>' 



648 



ln(3)vr 
+ 8964 



IT 



58 r - 

TpNf 



C A 

C F a 2 
648 



+ 0(a 3 



(6.2) 



We have checked that the Landau gauge part of this agrees with \17\ [T8] but the a dependence 
differs since as we have noted we have renormalized a in accordance with [13]. For the vector 
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the situation is effectively trivial due to the Slavnov- Taylor identity and so to the order we have 
computed 

C v (a,a) = 1 + 0(a 3 ) . 



For the tensor operator we have 



C T (a,a) 



1 + 



(3a - 3)^'(|) + (2 - 2a)vr 2 - 3a + 12 



Cpa 



(216a 2 - 432a + 504)(^'(|)) 2 + (576a - 288a 2 - 672)V>'(§)7r 



// 1 



+ (13824a - 432a 2 - 40176)^'(§) + (144 - 108a)V>"'(±' 
+ (62208a - 233280)s 2 (f ) + (466560 - 124416a)s 2 (§) 
+ (388800 - 103680a)s 3 (§) + (82944a - 311040)s 3 (§) 
+ (96a 2 + 96a - 160)vr 4 + (288a 2 - 9216a + 26784)vr 2 
+ 216a 2 - 1728a - 45063 + (1944a - 1944)S 

+ (5184a + 23328)C(3) + (432a - 1620) — ^ 

v3 

lnf3W 7r 3 
+ (19440 - 5184a) K ' + (1740 - AMa)-= 
v 3 v 3 



C F 



+ 



192V''(i)vr 2 - 144(V/(i)) 2 + (486a 2 - 1620a + 13590)^ 
+ (81 - 27a)V> w (!) + (124416 - 23328a)s 2 (f ) 



'(V 



+ (46656a - 248832)s 2 (f) + 



a - 207360)s 3 (f; 



+ (165888 - 31104a)s 3 (§) + (72a - 280)vr 4 

- (324a 2 - 1080a + 9060)vr 2 - 486a 2 - 2268a + 76419 



+ (972a - 1620)S - 32400C(3) + (864 - 162a) 



ln 2 (3)^ 



+ (1944a - 10368) + (174a - 928)-^= 



C A 



+ 



1440V- (|) - 960vr 2 - 17028 T F N f 



C F a 2 
1944 



+ 0(a 3 



where C T (a, 0) agrees with \17\ 118]. Numerically we have 



(6.3) 



(6.4) 



C s (a,a) = 1 + [0.2292715a - 0.6455188]a 

+ [0.5684263a 2 + 4.5546643a - 22.6076874 + 4.0135395iV / ]a 2 + 0(a 3 ) 

C v {a,a) = 1 + 0(a 3 ) 

C T (a,a) = 1 + [1.1181604a + 0.2151729]a 

+ [3.7661435a 2 + 10.8071579a + 43.4302495 - 4.1032786iV / ]a 2 

+ 0(a 3 ) . (6.5) 

For the tensor current in order to compare with the alternative scheme we have the alternative 
scheme conversion function to MS 



C T {a,a) 



alt 



1 + 

+ 



(15a - 15)V>'(|) + (10 - 10a)vr 2 - 27a + 54 



C F a 



1800a 2 - 3600a + 4392)(?//( ±)) 2 + (4800a - 2400a 2 - 5856)V>'(§)vr 2 
+ (58320a - 6480a 2 - 137376)^' (±) + (1296 - 540a)V>'"(±) 
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+ (186624a - 513216)s 2 (f ) + (1026432 - 373248a)s 2 (§) 
+ (855360 - 311040a)a 3 (f ) + (248832a - 684288)s 3 (f ) 
+ (800a 2 - 160a - 1504)vr 4 + (4320a 2 - 38880a + 91584)tt 2 
+ 5832a 2 - 11664a - 133893 + (9720a - 9720)£ 



+ (15552a + 38880)C(3) + (1296a - 3564)^-^ 

v3 

N ln(3)vr . . vr 3 

15552a) y + (3828 - 1392aW 
v3 v 3 



(42768 



C F 



+ 



1728-0'(±)7r 2 - 1296(V' / (3)) 2 + ( 24 30a 2 - 324a + 16542)^ 

- (135 + 135a)^'"(|) + (139968 - 69984a)s 2 (f ) 
+ (139968a - 279936)s 2 (§) + (116640a - 233280)s 3 (f ) 
+ (186624 - 93312a)s 3 (f ) + (360a - 216)vr 4 

- (1620a 2 - 216a + 11028)vr 2 - 4374a 2 - 20412a + 251397 

+ (4860a - 6804)S - 69984C(3) + (972 - 486a) - y* 

v3 

, ln(3)7T , , 7T 3 

+ (5832a - 11664) y + (522a - 1044)-= 
v3 v 3 



C A 



7200^ 



4800vr 2 



59724 



T F N f 



5832 



+ 0(a 3 



(6.6) 



Numerically this equates to 



C T (a,a) 



alt 



1 + [1.2710080a + 0.0623253]o 

+ [4.4752296a 2 + 12.2915908a + 42.4780444 - 3.9334478iV f ]a 2 

+ 0(a 3 ) . (6.7) 



In |19j it was noted that the Landau gauge one loop correction to the conversion function in this 
alternative scheme was significantly smaller than that of the scheme of [17[ 118]. However, with 
the two loop computations complete it transpires that the two loop terms are comparable. 

One of the roles of the conversion function is to allow one to map anomalous dimensions from 
one renormalization scheme to another. For instance, in the tensor case the two loop RI'/SMOM 
scheme anomalous dimension, (|5.6p . can be deduced from the MS version from C T (a,a) using 

TM'/SMOM (°RI' ' a RI') = TmS ( a Ms) ~ ^ ( a Ms) qIZZ ln ° T ( a MS' a Ms) 



MS 



~ «MS^ S ( a Ms) Q— ln ° T ( 



a MS' a MsJ 



(6.8) 



where in this case we have been careful in making it explicit what scheme the variables are in. 
We have checked that the RI'/SMOM scheme anomalous dimension of both the scalar and tensor 
follow precisely from the respective conversion functions. However, given the way the coupling 
constants appear throughout it is only the one loop part of the conversion function which is used 
in this. Therefore, given the fact that the three loop MS scalar and tensor operator anomalous 
dimensions are known then it is possible to deduce the three loop RI'/SMOM for the tensor 
current. We have verified that the three loop Landau gauge results of [T7] are obtained. 
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7 Discussion. 



We have computed the full two loop form of the Green's function of Figure 1 for three quark 
currents in both the MS and RI'/SMOM renormalization schemes. The amplitudes for the 
former case will assist with matching lattice computations of the same Green's functions at high 
energy once the numerical results are expressed in the same scheme. The latter amplitudes 
will play a similar role but in the case where the renormalization scheme used on the lattice is 
the same as the RI'/SMOM scheme defined in |161 [T71 118] . However, as we have discussed, in 
the case of the tensor current, there are in principle different ways of defining an RI'/SMOM 
type scheme. This is because with free Lorentz indices in the operator there is more than one 
amplitude with respect to a basis of Lorentz tensors. The tensor bases we have introduced here 
for the various operators are by no means unique. Since the renormalization constant is defined 
in relation to a basis tensor or tensors, then a change of basis would lead to different numerical 
structure in the renormalization constants themselves as well as the conversion functions. Whilst 
this degree of ambiguity in defining an RI'/SMOM scheme for these currents may appear to be 
an issue, it may in fact be possible to exploit it to render corrections in, say, the conversion 
functions such that they are significantly small. Therefore, one in principle can improve the 
accuracy of any measurements. For the tensor case it was noted in [19] that an alternative 
scheme in the tensor current case could produce a smaller one loop correction than that of 
the original work of [16|. [T71 118j . However, our two loop computation demonstrated that this 
is not retained at that order. Though one could conceive of a way of producing a smaller 
correction with the absorption of an appropriate finite part from another amplitude by a basis 
redefinition. Whilst we have focused on quark currents the programmes we have developed 
can be used to consider the same problem for the low moment operators used in deep inelastic 
scattering, [20]. The one loop analysis was given in [19] for the operator sets labelled W2 and 
W3. However, the treatment of the vector current here underpins any two loop extension of 
[19] . This is because within each of the sets W2 and W3 the vector current is present. Though 
as the deep inelastic scattering operators involve the covariant derivative the presence of the 
vector current is through its divergence since there is operator mixing. As a result of this the 
renormalization of the resident vector current still has to be consistent with the Slavnov- Taylor 
identity in the RI'/SMOM renormalization scheme. Therefore, we have laid the groundwork for 
that in this article by treating the vector current and its amplitudes when inserted in a quark 
2-point function at length. 
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Table 1. MS and RI'/SMOM scheme coefficients of Cf for one loop S amplitudes. 
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Table 2. MS and RI'/SMOM coefficients of CpTpNj for two loop S amplitudes. 
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Table 3. MS and RI'/SMOM coefficients of C F C A for two loop S amplitudes. 
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Table 4. MS and RI'/SMOM coefficients of C 2 F for two loop S amplitudes. 
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Table 5. MS coefficients of Cf for one loop V amplitudes. 
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Table 6. MS coefficients of CfTfNj for two loop V amplitudes. 
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Table 7. MS coefficients of CfCa for two loop V amplitudes. 
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Table 8. MS coefficients of C\ for two loop V amplitudes. 
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Table 9. RI'/SMOM coefficients of C F for one loop V amplitudes. 
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Table 10. RI'/SMOM coefficients of C F T F N f for two loop V amplitudes. 
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Table 11. RI'/SMOM coefficients of C F C A for two loop V amplitudes. 
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(V '(1/3)) 2 


4/9 


16/27 


8/27 


8/9 


/ /o\\2 

(V> (1/3)) a 














(V/(l/3)) 2 a 2 














^"'(1/3) 


1/6 


8/27 


1/27 


1/27 


^"(l/3)a 


-2/27 


-2/27 


-2/27 





vr 3 a/v% 
vr 3 /V3 


-29/243 





-58/243 





-29/486 


-29/81 


58/243 


29/27 


7rln(3)a/\/3 


-4/3 





-8/3 





vrln(3)/^/3 
vr(ln(3)) 2 a/^/3 


-2/3 


-4 


8/3 


12 


1/9 





2/9 





vr(ln(3)) 2 /V3 


1/18 


1/3 


-2/9 


-1 



Table 12. RI'/SMOM coefficients of C\ for two loop V amplitudes. 
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TT 2 a 
a 

IT 2 



^(1/3) 
y/(l/3)a 



TTT) 



MI 



TTT) 



TTT) 



( 1 ^ Tl TV Tfl 

« MS 


i z ; n MS 


W ra MS 


c 

W n MS 


MS 


W n MS 


2 





4/3 


2/3 


8/3 





-10/27 


-4/27 


-8/27 


-4/27 


-16/27 


4/27 


-2 





-4/3 


-2/3 


-8/3 





10/27 


-4/9 


8/27 


4/27 


16/27 


4/9 


-5/9 


2/3 


-4/9 


-2/9 


-8/9 


-2/3 


5/9 


2/9 


4/9 


2/9 


8/9 


-2/9 



Table 13. MS coefficients of Cf for one loop T amplitudes. 



21 


X(21) 
C (l)n 


MS 


^,(21) 
C (2)n 


MS 


C (3)n 


MS 


x(2i) 

C (4)n 


MS 


;j',(2i) 

C (7)n 


MS 


C (8)n 


MS 


1 


-182/27 







-80/27 


-40/27 


-160/27 







7r 2 a 
































a 
































7T 2 


-200/243 


128/81 


-160/243 


-80/243 


-320/243 


-128/81 


^(1/3) 


100/81 


-64/27 


80/81 


40/81 


160/81 


64/27 


Y>'(l/3)a 

































Table 14. MS coefficients of CpTpNf for two loop T amplitudes. 
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22 


T,(22) 




r,(22) 




T,(22) 




r,(22) 




'1 '.(22) 




'1 '.(22) 




ok ' 


C (l)n 


MS 


C (2)n 


MS 


C (3)n 


MS 


C (4)n 


MS 


C (7)n 


MS 


C (8)n 


MS 


1 


7097/216 







205/27 


205/54 


AW/27 







ir 2 a 


1/27 


16/27 


-28/27 


-14/27 


-56/27 


-16/27 


7r 4 a 


5/81 


2/27 


4/81 


2/81 


8/81 


-2/27 


C(3)a 




3 




2/3 



















2/3 


Sq 




5/6 




1/3 




2/3 




1/3 


4/3 




1/3 


a 




-10 







-14/3 




7/3 


-28/3 







vr 2 a 2 


-5/18 


-5/27 




2/9 




1/9 




4/9 


5/27 


a 2 


-15/8 









-1 




1/2 




-2 







vr 2 a 3 
































a 3 
































7T 2 


-919/486 


503/81 


1346/243 


673/243 


2692/243 


-503/81 


7T 4 


-1/27 




2/3 


52/243 


26/243 


104/243 




2/3 


C(3) 




-9 


26/3 


28/3 


14/3 


56/3 


-26/3 


s 




7/6 




5/3 




2/3 




1/3 




4/3 




5/3 


s 2 (7r/6) 




24 




112 




-80 




-40 




160 




112 


S2(7r/6)a 




-12 




-16 



















16 


S2 (vr/2) 




-48 




224 




160 




80 




320 




224 


s 2 (vr/2)a 




24 




32 



















-32 


s 3 (vr/6) 




-40 


560/3 


400/3 


200/3 


800/3 


-560/3 


S3(7r/6)a 




20 


80/3 

















-80/3 


S3 (vr/2) 




32 


-448/3 


-320/3 


-160/3 


-640/3 


448/3 


S 3 (7t/2)q! 




-16 


-64/3 

















64/3 


^(1/3) 


919/324 


-503/54 


-673/81 


-673/162 


-1346/81 


503/54 


^'(l/3)a 


-1/18 




8/9 


14/9 




7/9 


28/9 




8/9 


*//(l/3)ar 


5/12 


5/18 




1/3 




1/6 




2/3 


-5/18 


/ON 3 

?/) (l/ojar 
































^(l/3)vr 2 


8/27 


-16/27 


32/81 


16/81 


64/81 


16/27 


(V'(l/3)) 2 




2/9 




4/9 


-8/27 


-4/27 


-16/27 




4/9 


y/"(l/3) 


-5/216 


-19/108 


-7/54 


-7/108 


-7/27 


19/108 


^'"(l/3)a 


-5/216 


-1/36 


-1/54 


-1/108 


-1/27 


1/36 


7r 3 a/v% 
vr 3 /^3 


29/324 


29/243 

















-29/243 


-29/162 


203/243 


145/243 


145/486 


290/243 


-203/243 


7rln(3)a/\/3 




1 




4/3 



















4/3 


7rln(3)/V3 
vr(ln(3)) 2 a/^/3 




-2 


28/3 


20/3 


10/3 


40/3 


-28/3 


-1/12 




1/9 



















1/9 


vr(ln(3)) 2 /V3 




1/6 




7/9 




5/9 


-5/18 


-10/9 




7/9 



Table 15. MS coefficients of CfCa for two loop T amplitudes. 
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23 


r,(23) 




7 ',(23) 




r,(23) 




r,(23) 




r,(23) 




r,(23) 




an 


C (l)n 


MS 


C (2)n 


MS 


C (3)n 


MS 


C (4)n 


MS 


C (7)n 


MS 


C (8)n 


MS 


1 


-79/3 









-4 




-2 




-8 







-K 2 a 


-110/27 


-100/27 


-16/27 


-8/27 


-32/27 


100/27 


7r 4 a 


20/81 


8/81 


16/81 


8/81 


32/81 


-8/81 


C(3)a 




8/3 




8/3 



















8/3 


Sq 




5/3 




2/3 




4/3 




2/3 




8/3 




2/3 


a 




4 







16/3 




8/3 


32/3 







vr 2 a 2 


-10/27 







-8/27 


-4/27 


-16/27 







vr 4 a 2 
































a 2 




-1 









4/3 




2/3 




8/3 







7T 2 


128/9 


-148/9 


56/27 


28/27 


112/27 


148/9 


7T 4 


-32/81 


-16/81 


-128/243 


-64/243 


-256/243 


16/81 


C(3) 


20/3 




8/3 


-32/3 


-16/3 


-64/3 




8/3 


E 




5/3 




2 




4/3 




2/3 




8/3 




-2 


S2 (vr/6) 




-88 




128 




64 




32 




128 




128 


S2(7r/6)a 




32 




32 



















-32 


S2 (tt/2) 




176 




256 




128 




-64 




256 




256 


s 2 (vr/2)o! 




-64 




-64 



















64 


s 3 (vr/6) 


440/3 


-640/3 


-320/3 


-160/3 


-640/3 


640/3 


s 3 (7t/6)q! 


-160/3 


-160/3 

















160/3 


s 3 (vr/2) 


-352/3 


512/3 


256/3 


128/3 


512/3 


-512/3 


s 3 (7r/2)a 
^(1/3) 
^'(l/3)7r 2 a 


128/3 


128/3 

















-128/3 


-64/3 


74/3 


-28/9 


-14/9 


-56/9 


-74/3 
































Y>'(l/3)a 


55/9 


50/9 




8/9 


4/9 


16/9 


-50/9 


^'(l/3)7r 2 a 2 
V>'(l/3)a 2 


































5/9 







4/9 




2/9 




8/9 







^(l/3)vr 2 


-16/27 


32/27 


-64/81 


-32/81 


-128/81 


-32/27 


(1/3)) 2 


4/9 




8/9 


16/27 


8/27 


32/27 




8/9 


l„iJ(-\ ln\\2„ 

\W (1/3)) a 
(V>'(l/3)) 2 a 2 






























































y/"(i/3) 




2/9 


-2/27 


8/27 


4/27 


16/27 


2/27 


^'"(l/3)a 


-5/54 


-1/27 


-2/27 


-1/27 


-4/27 


1/27 


vr 3 a/V^ 

7T 3 /V3 


-58/243 


-58/243 

















58/243 


319/486 


-232/243 


-116/243 


-58/243 


-232/243 


232/243 


7rln(3)a/\/3 




8/3 




8/3 



















8/3 


7rln(3)/V3 
vr(ln(3)) 2 a/^/3 


22/3 


-32/3 


-16/3 




8/3 


-32/3 


32/3 




2/9 




2/9 



















2/9 


vr(ln(3)) 2 /V3 


-11/18 




8/9 




4/9 




2/9 




8/9 




8/9 



Table 16. MS coefficients of C\ for two loop T amplitudes. 
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1 


T (1 
C (Dn 


T (1 

C (2)n 


T (1 

C (3)n 


T 1 

* A 1 -) 

(4)n 


T (1) 
„ J A 1 -) 

°(7)n 


T 1 

i A 1 -) 

C (S)n 


1 


2/3 





4/3 


2/3 


8/3 





TT 2 a 


-4/27 


-4/27 


-8/27 


-4/27 


-16/27 


4/27 


a 


-2/3 





-4/3 


-2/3 


-8/3 





7T 2 


4/27 


-4/9 


8/27 


4/27 


16/27 


4/9 


^(1/3) 


-2/9 


2/3 


-4/9 


-2/9 


-8/9 


-2/3 


^'(l/3)a 


2/9 


2/9 


4/9 


2/9 


8/9 


-2/9 



Table 17. RI'/SMOM coefficients of C F for one loop T amplitudes. 



(21) 


T,{21) 
°(l)n 


C (2)n 


C (3)n 


T,{21) 
C (4)n 


T,{21) 
°(7)n 


X(21) 
C (S)n 


1 


-40/27 





-80/27 


-40/27 


-160/27 





ir 2 a 


80/243 


80/243 


160/243 


80/243 


320/243 


-80/243 


a 


40/27 





80/27 


40/27 


160/27 





7T 2 


-80/243 


128/81 


-160/243 


-80/243 


-320/243 


-128/81 


^(1/3) 


40/81 


-64/27 


80/81 


40/81 


160/81 


64/27 


y/(i/3)a 


-40/81 


-40/81 


-80/81 


-40/81 


-160/81 


40/81 



Table 18. RI'/SMOM coefficients of C F T F N f for two loop T amplitudes. 
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22 


T (22) 


T (22) 
C (2)n 


T (22) 
C (3)n 


T (22) 
C (4)n 


T (22) 
C (7)n 


T (22) 
C (8)n 


1 


205/54 





205/27 


205/54 


410/27 





7r 2 a 


-223/243 


47/243 


-446/243 


-223/243 


-892/243 


-47/243 


4 

tt a 


2/81 


2/27 


4/81 


2/81 


8/81 


-2/27 


C(3)q 





2/3 











-2/3 


Sq 


1/3 


1/3 


2/3 


1/3 


4/3 


-1/3 


a 


-223/54 





-223/27 


-223/54 


-446/27 





vrV 


-5/27 


-7/27 


-10/27 


-5/27 


-20/27 


7/27 


a 2 


-5/6 





-5/3 


-5/6 


-10/3 





7r 2 a 3 


-1/27 


-1/27 


-2/27 


-1/27 


-4/27 


1/27 


a 3 


-1/6 





-1/3 


-1/6 


-2/3 





TT 2 


673/243 


503/81 


1346/243 


673/243 


2692/243 


-503/81 


TT 4 


26/243 


2/3 


52/243 


26/243 


104/243 


-2/3 


C(3) 


14/3 


26/3 


28/3 


14/3 


56/3 


-26/3 


E 


-1/3 


5/3 


-2/3 


-1/3 


-4/3 


-5/3 


s 2 (vr/6) 


-40 


-112 


-80 


-40 


-160 


112 


s 2 (tt/6)q 





-16 











16 


S2 (tt/2) 


80 


224 


160 


80 


320 


-224 


s 2 (7t/2)q 





32 











-32 


s 3 (tt/6) 


200/3 


560/3 


400/3 


200/3 


800/3 


-560/3 


s 3 (7t/6)q 





80/3 











-80/3 


S3 (vr/2) 


-160/3 


-448/3 


-320/3 


-160/3 


-640/3 


448/3 


S3(7r/2)Q 





-64/3 











64/3 


y/(l/3) 


-673/162 


-503/54 


-673/81 


-673/162 


-1346/81 


503/54 


^'(l/3)a 


223/162 


-47/162 


223/81 


223/162 


446/819 


47/162 


0'(l/3)a 2 


5/18 


7/18 


5/9 


5/18 


10/9 


-7/18 


■0 (l/3)cr 


1/18 


1/18 


1/9 


1/18 


2/9 


-1/18 


^(l/3)vr 2 


16/81 


-16/27 


32/81 


16/81 


64/81 


16/27 


(V'(l/3)) 2 


-4/27 


4/9 


-8/27 


-4/27 


-16/27 


-4/9 


<'(l/3) 


-7/108 


-19/108 


-7/54 


-7/108 


-7/27 


19/108 


Y>"'(l/3)a 


-1/108 


-1/36 


-1/54 


-1/108 


-1/27 


1/36 


vr 3 a/v% 
vr 3 /^ 





29/243 











-29/243 


145/486 


203/243 


145/243 


145/486 


290/243 


-203/243 


7r ln(3)a/ \/3 





4/3 











-4/3 


7rln(3)/V3 
vr(ln(3)) 2 a/\/3 


10/3 


28/3 


20/3 


10/3 


40/3 


-28/3 





-1/9 











1/9 


^(ln(3)) 2 /V3 


-5/18 


-7/9 


-5/9 


-5/18 


-10/9 


7/9 



Table 19. RI'/SMOM coefficients of C F C A for two loop T amplitudes. 
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23 


T,(23) 


T,(23) 


T,(23) 


T,(23) 


T,(23) 


T,(23) 


a> n ' 




C (2)n 


C (3)n 


C (4)n 


C (7)n 


C (8)n 


1 


-10/9 





-20/9 


-10/9 


-40/9 





7r 2 a 


-80/81 


-268/81 


-160/81 


-80/81 


-320/81 


268/81 


4 

7r a 


8/243 


40/2431 


16/243 


8/243 


32/243 


-40/243 


C(3)q 





8/3 











-8/3 


Set 


2/3 


2/3 


4/3 


2/3 


8/3 


-2/3 


a 


8/9 





16/81 


8/9 


32/9 





vr 2 a 2 


16/81 


16/81 


32/81 


16/81 


64/81 


-16/81 


vr 4 a 2 


8/243 


8/243 


16/243 


8/243 


32/243 


-8/243 


a 2 


2/9 





4/9 


2/9 


8/9 





vr 2 


112/81 


-460/27 


224/81 


26/27 


448/81 


460/27 


7T 4 


-56/243 


-8/27 


-112/243 


-56/243 


-224/243 


8/27 


C(3) 


-16/3 


8/3 


-32/3 


-16/3 


-64/3 


-8/3 


£ 


-2/3 


2 


-4/3 


-2/3 


-8/3 


-2 


s 2 (vr/6) 


32 


128 


64 


32 


128 


-128 


s 2 (7r/6)a 





32 











-32 


S2 (tt/2) 


-64 


-256 


-128 


-64 


-256 


256 


s 2 (7r/2)a 





-64 











64 


s 3 (vr/6) 


-160/3 


-640/3 


-320/3 


-160/3 


-640/3 


640/3 


s 3 (7r/6)a 





-160/3 











160/3 


S3 (vr/2) 


128/3 


512/3 


256/3 


128/3 


512/3 


-512/3 


s 3 (7r/2)a 





128/3 











-128/3 


y/(l/3) 


-56/27 


230/9 


-112/27 


-56/27 


-224/27 


-230/9 


^'(l/3)7r 2 a 


16/81 


-16/81 


32/81 


16/81 


64/81 


16/81 


^(l/3)a 


40/27 


134/27 


80/27 


40/27 


160/27 


-134/27 


^'(l/3)vr 2 a 2 


-8/81 


-8/81 


-16/81 


-8/81 


-32/81 


8/81 


^'(l/3)a 2 


-8/27 


-8/27 


-16/27 


-8/27 


-32/27 


8/27 


^(1/3)7T 2 


-40/81 


40/27 


-80/81 


-40/81 


-160/81 


-40/27 


(y/(i/3)) 2 


10/27 


-10/9 


20/27 


10/27 


40/27 


10/9 


(V> (1/3)) a 
0//(l/3)) 2 a 2 


—4/27 


4/27 


—8/27 


—4/27 


—16/27 


—4/27 


2/27 


2/27 


4/27 


2/27 


8/27 


-2/27 


<'(l/3) 


4/27 


-2/27 


8/27 


4/27 


16/27 


2/27 


^"(l/3)a 


-1/27 


-1/27 


-2/27 


-1/27 


-4/27 


1/27 


7r 3 a/v% 





-58/243 











58/243 


-58/243 


-232/243 


-116/243 


-58/243 


-232/243 


232/243 


7r ln(3)a/\/3 





-8/3 











8/3 


7rln(3)/>/3 
vr(ln(3)) 2 a/^/3 


-8/3 


-32/3 


-16/3 


-8/3 


-32/3 


32/3 





2/9 











-2/9 


vr(ln(3)) 2 /V3 


2/9 


8/9 


4/9 


2/9 


8/9 


-8/9 



Table 20. RI'/SMOM coefficients of C 2 for two loop T amplitudes. 
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„ 23 


7 ',(23) 




r,(23) 




T,(23) 




'1 '.(23) 




7 ',(23) 




a n 


C (2)n 


alt 


C (3)n 


alt, 


C (4)n 


alt 


C (7)n 


alt 


C (8)n 


alt, 


1 







-4 


1/3 


r 


!/3 


-8/3 







it 2 a 


-28/9 


-224/81 


-112/81 


-448/81 


28/9 


it 4 a 


152/729 


-16/729 


-8/729 


-32/729 


-152/729 


C(3)a 


8/3 

















-8/3 


Sa 


2/3 


4/3 


2/3 


8/3 


-2/3 


a 



























vr 2 a 2 


8/27 


64/81 


32/81 


128/81 


-8/27 


vr 4 a 2 


40/729 


80/729 


40/729 


160/729 


-40/729 


a 2 







4/3 


2/3 


8/3 







7T 2 


-52/3 


256/81 


128/81 


512/81 


52/3 


7T 4 


-88/243 


-304/729 


-152/729 


-608/729 


88/243 


C(3) 


8/3 


-32/3 


-16/3 


-64/3 


-8/3 


s 




2 


-4/3 


-2/3 


-8/3 




-2 


s 2 (vr/6) 


128 




64 




32 


128 


-128 


S2(vr/6)a 




32 



















32 


S2 (vr/2) 


-256 


-128 




64 


-256 


256 


s 2 (7r/2)a 




64 



















64 


s 3 (vr/6) 


-640/3 


-320/3 


-160/3 


-640/3 


640/3 


s 3 (7r/6)a 


-160/3 

















160/3 


S3 (vr/2) 


512/3 


256/3 


128/3 


512/3 


-512/3 


s 3 (7r/2)a 


128/3 

















-128/3 


y/(l/3) 
^'(l/3)vr 2 a 




26 


-128/27 


-64/27 


-256/27 




26 


-80/243 


160/243 


80/243 


320/243 


80/243 




14/3 


112/27 


56/27 


224/27 


-14/3 


^'(l/3)7r 2 a 2 
^'(l/3)a 2 


-40/243 


-80/243 


-40/243 


-160/243 


40/243 


-4/9 


-32/27 


-16/27 


-64/27 


4/9 


^(1/3)tT 2 


136/81 


-272/243 


-136/243 


-544/243 


-136/81 


(V>'(l/3)) 2 


-34/27 


68/81 


34/81 


136/81 


34/27 


(ip (1/3)) a 


20/81 


-40/81 


-20/81 


-80/81 


-20/81 


(V>'(l/3)) 2 a 2 


10/81 


20/81 


10/81 


40/81 


-10/81 


<'(l/3) 


-2/27 


8/27 


4/27 


16/27 


2/27 


^"'(l/3)a 


-1/27 


-2/27 


-1/27 


-4/27 


1/27 


vr 3 /V3 


-58/243 

















58/243 


-232/243 


-116/243 


-58/243 


-232/243 


232/243 


7r ln(3)a/ \/3 


-8/3 

















8/3 


7rln(3)/>/3 
vr(ln(3)) 2 a/^/3 


-32/3 


-16/3 


-8/3 


-32/3 


32/3 


2/9 

















-2/9 


vr(ln(3)) 2 /V3 


8/9 


4/9 


2/9 


8/9 


-8/9 



Table 21. Alternative RI'/SMOM scheme coefficients of C\ for two loop T amplitudes. 
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